Entanglement and spin squeezing of Bose condensed atoms 



(N 
O 
O 
(N 

O 

Q 

(N 



> 
oo 

O 



o 



M. Zhang^, Kristian Helmerson^, and L. You^ 
^School of Physics, Georgia Institute of Technology, Atlanta, GA 30332-0430 and 
^"Atomic Physics Division, National Institute of Standards and Technology, Gaithershurg, MD 20899-8424 

(Dated: February 1, 2008) 

We propose a direct, coherent coupling scheme that can create massively entangled states of 
Bose-Einstein condensed atoms. Our idea is based on an effective interaction between two atoms 
from coherent Raman processes through a (two atom) molecular intermediate state. We compare 
our scheme with other recent proposals for generation of massive entanglement of Bose condensed 
atoms. We also present explicit calculations that support detecting maximally entangled states 
without requiring single atom resolved measurements. 



I. INTRODUCTION 

Entanglement lies at the heart of the difference be- 
tween the quantum and classical multi-particle world, ft 
is the phenomenon that enables quantum information 
and quantum computing with many qubits. Recently, 
several interesting developments have occurred in stud- 
ies of massively entangled atomic states. Based on the 
proposals of Sorensen and Molmer a controlled, en- 
tangled state of 4— ions was successfully created by the 
NIST ion trap group Zeilinger and coworkers pre- 
pared three entangled photon or GHZ states by selecting 
from two beams of entangled photon pairs |^ . Entangle- 
ment between two atoms and a microwave photon were 
also detected in a "step-by-step" process 

Of these and other related developments, the idea of 
Molmer and Sorensen [Q is especially interesting. In their 
scheme a direct coupling to the multiparticle, entangled 
final state was made possible through a virtual, interme- 
diate state which was a common (quantum) mode of the 
motion of all the ions. Similar type interactions were also 
obtained by Milburn Both proposals allow for cre- 
ation of massive entangled states by unitary evolution, 
starting from certain pure initial states. 

A Bose-Einstein condensate is a convenient source 
of atoms in initially pure and separable states. ZoUer 
and coworkers suggested creating massive entangled spin 
squeezed states from a two component condensate us- 
ing the inherent atom-atom interactions |^ and further 
investigated schemes to dynamically create many parti- 
cle entangled states of a two component BEG 0. Spin 
exchange interaction of a spinor condensate was also pro- 
posed as a candidate for creating entangled pairs of atoms 
iJ. 

Most of these proposals work in the two mode ap- 
proximation where one motional state is assumed for 
each spinor component of condensed atoms. Recently, 
Sorensen demonstrated the validity of this method within 
Bogoliubov theory pO[ |. Similar studies were previously 
performed, extensively, for condensate atoms in a double 
well potential |ri|, . 

Recently, two of us proposed a new type of atom- 
atom coupling that achieves optimal spin squeezing in 
two mode Bose condensed atoms In this study 



we provide further details of our proposal and compare 
its potential for both the creation of massive entangled 
states Q and improved spin squeezing ||, ^, 18|. This 
paper is organized as follows. First, we briefly review 
the Raman process with photon momentum transfer for 
a single three-level A-type atom. We then extend the 
analysis to the case of two atoms, using intermediate 
molecular states. The resulting effective interaction is 
then compared with models in several other recent stud- 
ies of creating entanglement with Bose condensed atoms. 
We then provide an explicit calculation of the time re- 
versed dynamics to confirm the generation of maximally 
entangled states without requiring single atom resolved 
measurements |l^ . We conclude with a careful analy- 
sis of the limitations and advantages of our model. 



II. SINGLE ATOM, TWO-PHOTON PROCESS 
THROUGH AN ATOMIC INTERMEDIATE 
STATE 

We consider a three level A-type atom described by the 
Hamiltonian 



where the dipole interaction term is 



|Ai)(Ai|+Hi„t, (1) 



Hint = -~ig-Ei{r)e'^''''~''^''\9){e\ + h.c. 

~d,g, ■ J?2(r)e^"~-^~-^"=* \g') (e| + h.c. (2) 

for the general Raman process from initial state \g) to 
final state \g') through an intermediate state |e). V^(r) 
denotes the, possibly different, trapping potentials for 
the different electronic states with internal energy 
Ej{f)e'^^^'^~'^'^^* denotes the running 
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wave amplitude of the laser field. 

We define A = LOi—LOeg and 5j = 
of generality, assuming deg = deg' 



— ijji. Without loss 
d, nj 2d - Ej/h 



is the Rabi frequency for the corresponding dipole con- 
nected transition. In the limit of a resonant two photon 
process with large detuning from the intermediate atomic 
excited state |e), we adiabatically eliminate |e) to obtain 
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an effective Hamiltonian 



\tJ'){^^\ + VB/R, 



(3) 



with, 

Vb - 



4A 



wliere distinguish the two types Raman process according 
to whether the final internal state \g') is the same as or 
different from the initial state \g), respectively. When 
Id') — \9)^ identify the general Raman process as 
Bragg diffraction |Q . 

In momentum space, the above two photon process in- 
volves the simultaneous absorption of a photon (jti, uJij 
by an atom in state \g',p) and stimulated transition 
to state 



g';p — hk2 + tikij with emission of a photon 

(jt2,i^2^- For the case of Bragg diffraction {g' — g) 
with two counter-propagating waves 



fci = k, 
k2 « -fc, 
K = ki-k2 



2k, 



(6) 



we obtain 

VB{r,t) 



h 

4A 



iSot AK r 



\9){9\- (7) 



The first two terms are AC Stark shifts that can be ne- 
glected. For an atom initially in state \g,pg), the last two 

terms of Eq. (|^) couples to momentum states pg ± 2hk, 
respectively. By appropriately choosing 82 one can selec- 
tively enhance coupling to only one momentum state or 
side mode. The resonance condition is defined by energy 
conservation in the two photon process 



huji 

which gives 
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hK 



2M 



(8) 
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(9) 



for an atom with pg 0. Hence under conditions reso- 
nant for Bragg diffraction, the effective coupling (Q) cre- 
ates superposition states like. 



\i;)^a\g)\pg)+(3\g) pg±2hky (10) 

For Bose-condensed atoms, the above discussion still 
applies to any of the identical single atoms. When a 
condensate of N atoms is involved (initial motional state 
\pg ~ 0)), one simply creates a condensate in state (pX|), 
which is not an entangled state because it is simply 
putting N atoms into the same state (|l0|), i.e. 

(af^J |vac) ^ ■ --Wn) - IV^)^^ . (11) 

In practice this can only be done if atoms are nonin- 
teracting, or that the effect of interactions is small dur- 
ing the time of atomic Bragg diffraction. Bragg diffrac- 
tion of condensate atoms was demonstrated first by the 
NIST-Gaithersburg group |Q. Higher motional states 
K — ±2mk (to integer) was also obtained through, 
higher order, 2TO-photon couplings. 

For the Raman process with two co-propagating 
waves fci « k2 involving nearly degenerate ground states 
\g) and \g'), we have 



Vr 



h 

4A 



n,n*2e^'^' \g){g'\ 



(12) 



where AC Stark shifts from second order processes in- 
volving the same laser fields have been neglected. The 
Raman resonance condition is then simply wi— 0^2 = ^g'g- 
Note that co-propagating Raman results in an effective 
coupling (|i^) which is constant, independent of the posi- 
tion f. The two photon process then is essentially insen- 
sitive to motional effects {k2 ~ fci). 



III. TWO-ATOM, TWO-PHOTON PROCESS 
THROUGH AN INTERMEDIATE MOLECULAR 
STATE 



We now consider a model involving two A-type atoms 
whose initial and final states are described by the same 
non-interacting atomic states. The intermediate excited 
state, however, is now chosen to be a bound, molecular, 
excited state, similar to those utilized in recent photo- 
association experiments |2^. Neglecting configurations 
not directly involved in the two photon process, the two 
atom Hamiltonian can be written as 
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n = 



E 



1=1,2 



E (vb + ^s'airr)) + Ug^giln -r2\) + nug,g 

1=1,2 ^ ^ 



l5',5>(5',5l 



E + (^1)) + C^e. - I) + 

1=1,2 ^ ^ 



|e6)(e6| + Hint, 



(13) 



where states ^) for /i = g, g' denote symmetrized elec- 
tronic states of two atoms, each in \fj,). \g',g) denote 
the symmetrized state with one in \g) and the other in 
\g'). Vf^i^g^gi ^gi g^et dcnotc Bxtemal trapping potentials in 
state l/i), J7^^ and Ugig are atom-atom interactions which 
at large inter-atomic distances are described by the usual 
van der Waals type terms, and at short range modified by 
coulomb effects. The intermediate molecular state fam- 
ily, denoted by |ef,), can be accessed through the direct 
dipole coupling term iJint- C^e(, is the Born-Oppenheimer 
molecular (of atom 1 and 2) binding potential for inter- 
nal state manifold |ef,), which contain bound levels to be 
used as intermediate states. 

Introducing center of mass R — (ri +r2)/2 and relative 
coordinate r = ^2 — ri for the two atoms, we can express 



^ r ^ ^ r 
ri=R-^, r2 = R+-, 



and 



Pi , P2 



2M 2M 2(2M) 2(Af/2)' 



(14) 



(15) 



When \eb) is asymptotically connected with |e,/z) = 
(|e) l/z) -I- III) |e)) /-\/2 (for /i = g,g'), we can express the 
dipole coupling as 

iJint = - E • E^e^'^'-'^^-'^'* (|g)(e|), + h.c. 

i=l,2 

- 5^ • ^26*^"-^"-'-^* (|g')(e|). + h.c. 

1=1,2 



"7f 



cos(fci • f/2)e*'=i-^-*"i* \g,g){eb\ + h.c. 
cos(fc2 • r72)e*^-=-^-*"=* \g\g'){eb\ + h.c, (16) 



where we have assumed di = d'^ = d for the \e) 
transition of atom i. We note that the electronic dipole 
coupling f2j will vary considerably with r for small values 
of r, due to the relatively short range of the molecular 
interactions. 

To express the intermediate state in terms of the par- 
ticular molecular resonance or bound state le;,), we now 



take a look at the eigen-structure within the jcf,) fam- 
ily. We will essentially base our discussion on some kind 
of molecular state asymptotically connected to the limit 
of one ground \g) and one excited atom |e), just as the 
Og state, extensively discussed by D. Heinzen and others 
in their photo association work The recent exper- 

iment by Heinzen and coworkers on the production of 
ground state molecules from an atomic condensate by a 
two-photon Raman process pl| ] provides additional mo- 
tivation to explore our ideas experimentally. The photo- 
association process [^2[ |3[ ^J] used by Heinzen et al. 
relies on the transition strength of going from a (two 
atom) 'free' to a (molecular) 'bound' state via an inter- 
mediate, excited (molecular) 'bound' state. What we de- 
sire, on the other hand, is a transition from a (two atom) 
'free' to 'free' (relative motion) state via an intermediate, 
excited (molecular) 'bound' state. More detailed discus- 
sions about feasibility of such transitions can be given 
provided all molecular potential curves are available. 

The question of whether this is possible or not does not 
seem to depend on the sample density (in the weakly in- 
teracting limit), but will depend, as we shall show below, 
to a large degree on the trap strength and the excited 
bound state structure. Assuming all trapping potentials 
Vg, Vg> , and to be harmonic, the separation of center- 
of-mass and relative coordinates for the trapping poten- 
tials is 



J = l,2,3 j = l,2 

= I (2M) ^x] + i (M/2) 



2 2 
J = l,2.3 



(17) 



where 



Xij — Cj ■ R, and Xj = 



' with 



the unit vector in the j direction. 

When short range interactions U^^if) are approx- 
imated by their optical contact forms, Wilkens and 
coworkers provide analytic solutions for two interacting 
atoms inside a harmonic trap pq|. In principle, from 



P 



2(M/2) 



cj^ilif) ^ EXir), (18) 
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we can find all bound states 4>n{r) due to the external 
trapping potential U^^{r). We note that if the range of 
the inter-atomic interaction ^ r <C size of the ground 
state of the harmonic trap, then these bound states 
should resemble the standard low energy scattering so- 
lutions. 

Similarly, the bound state of the intermediate, excited 
molecular level satisfies 



2{M/2) 



-HV;?(f) + C/e,(f) 



(19) 



where the energy Ef^ is measured from the asymptotic 
energy hujeg- The index m represents the rotational and 
vibrational quantum numbers. For all bound states with 



J 



Ef-^ < 0, the effect of external trapping potential V^'^!' is 
much less than for the free atoms, since the corresponding 
bound state wave function is localized inside the molec- 
ular potential well Ue^ ■ 

Similarly, complete bound states |^) can also be intro- 
duced for \g',g) if needed. They satisfy 



2(M/2) 



Mr)^Ef'iji{r)- (20) 



In the {R, f) basis, keeping only one intermediate 
bound state mt, we can then write our previous Hamil- 
tonian, Eq. (|l3|), as 



+E 



p2 



2{2M) 



Vtj,{R) + 2hu^, + E^ 



^^+V^iR) + n.,, + Ef^ 



^^+Vr,(R) + n.^, + E:^^ 



9',9;i){g',9;i\ 

eb;mb){eb;mb\ + H-.^t, 



and the interaction term Eq. (16) now becomes 



E 



rigmi,'- 



' \9^9\ng){eb]mb\ + h.c. 



V2 

■"9 

E Vng,m,e'''-^^'~'^'' Iff', 5'; n,'> {eb; mb\ + h.c. 



r 



(21) 



(22) 



where 



Vn^m^ = -j{ri^\d{r)cos{k-r/2)\m.b) 
\d\ 



drcj)"^; {r)d{r) cos(fc • f/2)C.^(r), (23) 



which is a measure of the free-bound transition ampli- 
tude, we note that in contrast to other ex pre ssions for 
the free-bound transition amplitude ||2^, the mo- 

tional wave functions \n^) and \mb) are now normalized 
in the trap rather than energy or momentum normal- 
ized. Although the modification of the inner part of the 



J 



trap potential for r < 100 (oq) due to the inter-atomic 
potentials is to be included, these values can be deter- 
mined for selected molecular states from the results of 
photo-association experiments. They can also be com- 
puted directly if accurate potentials are available. De- 
tailed discussions are given by Verhaar |2^ , Julienne |2^ 
and Javaninen p^ . 

In general, the coupling coefficient will vary depending 
on intermediate state. It is therefore important to pick 
the largest one of them. We then perform an adiabatic 
elimination of the excited state by appropriate choice of 
laser detunings, and obtain the effective Hamiltonian, 



^-EE 
E 



2(2M) 



p2 



2(2Af) 



V,%%R) + nu,,, + E^, 



9',9'J){9',9;1\ + Vb/r, 



(24) 
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with 



4A 



^ikj -R—iujjt 



\g,g;n){g,g;l\ , 



ik2- R—iuJ2t 



W ■,9'\n'){g,g;n\ + h.c. 



(25) 



r 



where we have defined A = uji — L02 — E'^^/fi. Note the 
(v^)^ enhancement due to excitation to an asymptotic, 
symmetric state. 

More generally, one can include all motional states and 
still be able to find the effective coupling between the 
selected electronic states, to obtain similar results |26 . 



Our analysis now follows along the same lines as in the 
previous case of a single atom. For example, if we con- 
sider only the lowest motional state = (a situation 
well approximated by Bose condensed atoms), then our 
Hamiltonian becomes 



+ Vb/r, 



(26) 



with 



Ak-j -R—iijjj t 



R~iujj t 



|g,g;0)(ff,5;0'|, 



,ik2-R—iuJ2t 



t 



\g',g';0){g,g;0'\ + h.c. 
For Bragg diffraction by two counter- propagating waves, the above result reduces to 

'e~'^-^] |5,5;0)(5,ff;0| 
I 



VBiR,t) = ^ \r,o^S + 1^2^ + ni^le^'^'e'" '' + ^l^l^e 



(27) 



(28) 



Similar to the single atom case, Bragg diffraction involves 
the simultaneous absorption and stimulated emission of 
photons, i.e. the elementary Bragg process involves a to- 
tal of two photons. In contrast, however, to the single 
atom case, a pair of atoms are now diffracted as illus- 
trated in Fig. |l|. Hence for Bose condensed atoms with 
Pi ~ 0, Bragg diffraction produces atoms in the state 



a |g, g- 0) bi « 0, P2 « 0) + /3 I5, g\ 0) e±*^-« 
a \9.9) \Pi ~ 0,p2 ~ 0) 



+li\g.g) Pi = nK/2,p2 = nK/2) 



(29) 



Not surprisingly, the momentum shift is only half the 
value of atomic Bragg diffraction, due to the fact that 
only two photons are involved for the pair of atoms. The 
resonance condition occurs at 



huji 



2(2M) 2(M/2) 



p2 



+ 



2(2M) 2(M/2)^ 



Pi Pf ^0, 

Pf = P,±hK, 

which gives the resonance condition 



ns2 



{P^ ± TlKf 

2{2M) 2{2M) 



'2(2M)' 



(30) 



(31) 



i.e. half the atomic Bragg resonance (||). 

It is worth mentioning that for a condensate with many 
atoms, the state obtained from diffraction via a two atom 
intermediate state will be complicated. It is generally not 
of the form (aj^^ )^/^ |vac). 

For a nearly degenerate Raman process with two co- 
propagating waves, we obtain for atoms initially at rest. 



VR{R,t) 



2 ^ 'VOmtVo'rrib 

e''^'\g,g;0){g',g';0'\ + h.c. 



(32) 
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molecular excited state 
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FIG. 1; The illustration of energy level diagrams for Bragg 
diffractions through atomic (left part) and molecular excited 
states (right part). 



The resonance condition now becomes toi — uj2 — 'icogig, 
i.e. the Raman process changes each atom's state (g — s- 
g') and so each atom acquires the energy deficit Tiujgtg. 

Note that there is no R dependence of the coupHng coef- 
ficient in this case. 

If the Raman process is arranged to allow for simulta- 
neous scattering to two diff'erent final states \g') and 
the coherent coupling becomes, 



,1(52* 



2 



\g,g-0){g',g'-0'\ + h.c. 



\g,g-Q){g",g"-Q"\ + h.c. (33) 



Such a coupling can be obtained among different Zeeman 
levels in the atomic ground state. 

Finally, we note that our discussions above can also be 
applied to a pair of different species of atoms, e.g. using 
molecular intermediate state from Li-Cs dimer could cre- 
ate entangled pairs of Li and Cs. We also note that to be 
formally correct, the second order perturbation process 
also needs to be considered for the state which 
should have some type of effective couplings with itself 
and as well as with the states \g,g) and \g' ,g')- we have 
neglected those couplings, assuming that resonance con- 
ditions for their survival in the interaction picture are not 
satisfied. 



IV. SPIN-SPIN INTERACTION MODELS AND 
TWO MODE APPROXIMATIONS 



In this section, we will be primarily concerned with 
process which produce entanglement of atomic spin 
states. We therefore consider a spin 1/2 system and 
adopt the notation |t) for spin up and \[) for spin down, 
with the following operators and identities: 

cr. - ITXTl-liXil, 

= 2^^x+i(ry) = |t)(i| , 

^{<^X -iCTy) = U)(T|, 



1, 



(34) 



The Bragg/Raman coupling as discussed here is of the 
form 



[(IT)ai)i®(IT)ai)2 + (U)(TI)i®(U)(TI)2] 



^^R, (1) 



J2) 



+ (7i'> (35) 



(2)^ 



if we designate \g) and \g') as |t) and |J,), respectively. 
The Rabi frequency flu is the two-photon Rabi frequency 
for transitions between jt) and |J,) via the intermediate 
molecular state. Hence it contains the single photon Rabi 
frequencies fii and Q21 a-s well as the free-bound transi- 
tion amplitude, rjn^mf, ■ This is a new type of coupling not 
widely discussed before. Its ability, however, to generate 
entanglement should be obvious as two atoms perform 
conditional evolution at all times, similar to the coupling 
for photon down conversion. Its prospects for creating 
massive GHZ states will now be studied. 

We first compare our coupling scheme with other rel- 
evant models. In the original scheme of Sorensen and 
Molmer [Q, the two atom coupling takes the following 
form. 



nn 



r(2) 



(IT>ai + U>(TI)i®(IT> ai + U>(TI)2,(36) 



a form different from ours. It is convenient to analyze 
such models for the case of many atoms in terms of col- 
lective spin operator Jfj,=x,y,z = ^ X)i "'m''- One can then 
show that the Hamiltonian for the Sorensen and Molmer 
scheme becomes 



E 

i<j 



while our model gives 



+ 



N 



(37) 



.0) 



^^R i 72 72 ^ 



{Jx Jy) ■ 



(38) 



Recently, a many body, two mode coupling was pro- 
posed by Zoller and coworker |6|. They considered an 
interacting, two component (i.e. | [) and | ])) conden- 
sate whose many body Hamiltonian becomes 



H = 



7 



(39) 



under the two mode approximation, where each compo- 
nent has the same spatial mode, 



Introducing the coUective spin operators 



(40) 



Jx — 


1 

2 




+ ajao 


Jy — 


i 
2 




- a|ao 


J. = 


1 
2 




- ajao 



(41) 



and taking the index for | |) and 1 for | t) |T^, 
p7| , one can show that Jx,y,z = ^Ei^2;'?a,2 
used earher. We also note that N = ajai + a\aQ and the 
Casimir relation 



The interaction term Eq. ( 39 ) then becomes 
i J| (t/oo + C/ii - 2C/oi) /" dr |(/.o(f)|" , 



(42) 



(43) 



In the discussions of the numerical solutions to follow, 
we will use the second quantized notation, and consider 
a pure initial state with a fixed, total number of atoms 
iV. Specifically, we compare the three types of coupling. 



Vm = n^iRjl 



(44) 



considered by Sorensen and Molmer Q , the spin squeez- 
ing model 



Vs = nnnJl 



(45) 



recently considered by ZoUer and coworkers and our 
proposed coupling 



V, 



(46) 



Given an initial state |'!/'(0)) = |vac) /Vn\, i.e. 
with all atoms condensed in state | 1), the time evolu- 
tion operator for Vai is analytically known and takes a 
simple form at flfjt = 7r/2 for iV an even integer [Q, 

Um = e^^'^'^"^*) = e"^/2+"/4-latJV^^ ^ (47) 
It creates the maximally entangled N-GHZ state 

|GHZ)^ (X -L[|0)JiV),+.^+i|iV)jO)J. (48) 



In the context of Bose-Einstein condensation, state 
|GHZ)^ is an example of the interesting fragmented con- 
densate The Sorensen and Molmer model (Eq. (^)) 
produces perfect GHZ-type states at selected times, while 
our coupling (Eq. (^6|)) doesn't, in general, produce ex- 
act GHZ-type state. From numerical simulations, how- 
ever, we find that our model, on average, produces more 
than 50% overlap with GHZ-type state (Eq. (^)) at se- 
lected times. Because of the geometric equivalence be- 
tween and J|, the spin squeezing model (Eq. (|4^)) 
also produces a perfect GHZ at selected times, except 
now the entanglement is between single particle states 
|±) = (I t) ± I i))/\/2 |§ and the initial state has to 
become |±)®^. 

For our model, we were not able to find a general ana- 
lytic form of the time evolution operator, even at specific 
values of t. To numerically calculate the time evolution, 
we expand the wave function as 



N 



m)) = E^™wi^-"^)oi"^)i' 



where 



m=0 



m 



vac) 



(49) 



(50) 



with the initial conditions given by the Cm (0)s. 

In Figs. H and || we show, for comparison, the gener- 
ation of the maximum entangled GHZ state Eq. (Eq) as 
determined by model Vm (Eq. (0)) and our "molecu- 
lar diffraction" model (Eq. (p^)). We see that while the 
Sorensen and Molmer model creates perfect massive GHZ 
states at selected times, our model can also creates sig- 
nificant overlap with the massive GHZ at selected times. 

We can also compare the achievable spin squeezing be- 
tween our model and that of Vs (Eq. (^), using the 
squeezing parameter 



iV(AJsj2 



(51) 



where n^, i — 1,2,3 are mutually orthogonal unit vectors. 
Other discussions of spin squeezed states can be found in 
@|l|,||. 
For the initial state 



\m) = 



1 



1 



2^/2 



(al + a\Y \Ya.c), (52) 



it has been shown Q that $^^{t > 0) < 1 for some set of 
n^'s. This model of Vs (Eq. (|45|)) is in fact the one-axis 
twisting model considered by Kitagawa and Ueda earlier 
jlSf . In this case the problem can be solved analytically, 
with the result 



2 

mill 



(53) 



On the other hand, our model resembles the two-axis 
twisting model of Kitagawa and Ueda fl^, and has to 
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FIG. 2: The time dependent probabilities for being in states 
!0>ol^>i (|cjv|^) and |iV)o |0>i (|co|^) from the Vm model. The 
initial conditions are Cm(0) = 5mN {N — 500). For clarity, 
|cjvp (dashed line) is superimposed into the plot of |coP + 
Icjvj^. Note the perfect |GHZ)^ state at ^Ijit = 7r/2. 
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FIG. 4: Numerical results from the model Vs- (a) The time 
dependence of the minimal spin squeezing parameter; (b) The 
optimal angle when the minimal squeezing in (a) is achieved. 
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FIG. 3: The same as in Fig. ^ but for our model Vb/r- Note 
that at several places the projection onto the |GHZ)^ state is 
already significant. At fiflt ~ 0.22 (vertical line) the system 
exists in |GHZ)^ with about 50% probability. 




° 0.005 0.01 0.015 0.02 0.025 

FIG. 5: Same as in Fig. ^ but from our model Vb/r- Note 
that the minimal squeezing of our model is smaller and the it 
occurs at a constant spin direction. 



be solved numerically. In the limit of large N and with 
the condensate initially in | l)"*^^ one can show that for 
an initial Jz — N/2 (J-spin pointing along the positive 
z-direction), 

(AJsJ^i„^l (54) 

The optimal squeezing in this case occurs with Jx + Jy- 
This result can be easily verified by taking a semiclassi- 
cal approximation in the dynamical equation for and 
Jy. We find that the time scale of reaching maximum 
squeezing is ~ l/{Nflj^) (see also Q). For condensates 
containing 10^ atoms, even with a very weak coupling 



ilfi ~ 1 (Hz), the maximum squeezing is reached within 
a micro-second. 

In Figure ^ the squeezing parameters and the corre- 
sponding angle ^min SLie computed numerically for the 
model Vs (Eq. (|4^)) with the initial state given by 
Eq. ( |5^) . Fig. H is a similar calculation for our model 
(Eq. (|4(]|)) starting with all atoms in state | t)- Out 
model achieves a better squeezing at an earlier time. In 
addition, the phase angle ^min is fixed at the value 7r/4. 
More detailed discussions of the spin squeezing proper- 
ties of our model can beperformed, including dissipation 
and finite system size [|4[ |5[ |6|. 

Finally we also note that our coupling scheme for si- 
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multaneous transitions to two different final states gives 

r[2] T^^R 



B/R 



h.c. 



(55) 



In the undepleted pump limit when dissipation exists, or 
when the coupling efficiency is small, this coupling can 
be approximated by 



V, 



B/R 



£ ( a\a\ + a\,a\, 



(56) 



This is similar to the recent ideas of Refs. [|[ ^, ^ which 
creates entangled EPR pairs every time an event occurs, 
just like in a parametric down conversion process B^, B8[. 



V. DETECTION OF THE MAXIMALLY 
ENTANGLED STATES 

While it seems that there are promising schemes to 
generating massive entanglement of condensate atoms, 
the detection of such entanglement represents a signifi- 
cant challenge by itself as it is very difficult to perform 
individual atom resolved measurements in a con- 

densate. Recently, it was suggested that performing a 
time reversed dynamic evolution can detect the coher- 
ence of the created N-GHZ state 0, 

In this section, we detail the working mechanism of the 
proposed scheme by performing explicit calculations. We 
will use the Vm model as an example. To simplify the 
algebra, we will use r — ^Rt, and denote the basis states 
according to 

\ni) EE |no,ni) = |no)o l"-i>i • (57) 

Assume we start with | f)®^, then at r = 7r/2, a maxi- 
mally entangled state 

|V'(r = ^))(xi=(|7V)+z~+i|0)), (58) 

Immediately reversing the time dynamics (Vm —Vm), 
it is easy to show that after another t = 7r/2, we recover 
the initial state | l)"*^, i.e. all atoms to be detected are 
in state | t)- If instead, the created state at r = 7r/2 is 
a completely incoherent mixture of |0) and |iV), given by 
the (N-body) density matrix 



1 



1, 



|0)(0| + -|iV)(iV| 



(59) 



2' 2'"""' ' 2' 

then the time reversed dynamics will lead to observations 
of variable number of atoms in state | |) as well. 

For the general case, let's assume the restricted form 

p(r = |) = po\Om+PNmN\ 

+m{N\+l3*\Nm, (60) 

with < po/Af < ^, Po + Pn — 1, and < PqPn- Time 
reversed evolution of Vm for r from r = 7r/2 will lead to 



Pir + n/2) 



UM{-r)p{7r/2)Ul{~T) 
e^"''V(7r/2)e-*'^'^ 



(61) 



which can be explicitly evaluated with the use of rotation 
operation 



and the disentangling relation |39| 



Denote 



A ^ {N\e 
B = {N\ 



|0), 



(62) 



(63) 
(64) 



(65) 
(66) 



then we find 



Pnn{t + tt/2) = {N\p{t + tt/2)\N) 

= ^M^Po + \BfpN + {AB*I3 + h.c.)] (67) 



with 



m 



m\{N — to)! 



N 

E 

N 

^ m!(iV-TO)!' 

m=0 ^ ' 



(_-^)mg»(m-f)V^ (68) 



iV! 



(69) 



The following angular momentum algebra has been used. 



J+ |m) = ^(TO + l)(iV-m)|TO+l), 

J_ |to) = y^m{N - TO + 1) |to - 1) , 

Jz \m) : 

(to|J+ = (J_|m))t. (70) 



TO — I |TO) 



At r = f , we find 



PNN{7T) = - + {~1)-Ini{^)<l, 



(71) 



because |Im(/3)| < < y/poPN < 1/2, for any state. 
Thus we conclude that the probability of detecting atoms 
in state | J,) is nonzero for mixed states. In fact, for the 
completely incoherent mixture Eq. (|5^) with Po = Pn = 
\ and /? = 0, pmiir) = PArAr(7r) = 1/2. 

It's interesting to compare the average populations and 
fluctuations of spin up | ]) number operator [rii — a\ai) 
after the time reversed evolution for t — ^ starting with 
the following different initial conditions, 

• A pure state, 

|V'(7r/2)) = -i5(|iV)+z^+i|0)) 



(ni(^)) = tr(pni) = iV, 
(n2(7r)) = tv{pn\) = N\ 



= 0; 



(72) 
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FIG. 6: The time dependence of average population (ni = 
a\ai) and its variance (Ani) for spin up | |) component. The 
sohd hnes are for starting in the N-GHZ state Eq. ( |58[ ) while 
the dashed lines are for starting in the incoherent mixture Eq. 
(p9[). The upper two lines are for (ni) while the lower ones 
are for Ani. 

• An incoherent mixture, 
p(7r/2) = i(|7V) (iV| + |0) (0|) 

N 

(ni(7r)) = tr(pni) = —, 

{nl{n)) = tripnl) = — , 

Ani(7r) = (73) 

Thus, these two cases become easily distinguishable. In 
Fig. ^, we compare the time dependent results of aver- 
aged population in state | t) and its variance, from the 
time reversed evolution starting with the coherent GHZ 
state Eq (^) and the incoherent mixture Eq. (^). It 
is interesting to note that the average population of spin 
up I 1) number operator (ni — a\ai) for the completely 
incoherent mixture is in fact independent of r. We find 
that the width of the Gaussian like features near Qfit tt 
is inversely proportional to /Vn, but the plateau for the 
variance Ani is at around N/ ^/S. 

For factorizable N-atom states, whether a pure state 
IV-WS)) = 2W^(|T) + U))^"^ or a mixture p(7r/2) = 
27r(|i) (il + IT) (Tl)®^, the probability distribution after 
the time reversed evolution gives identical results that 
are different from the two multi-atom correlated states 
discussed earlier. In fact, the final results simply reflects 
the binomial distributions of rii in a equal superposition 
separable state and is given by 

N 

(ni{n)) = tr(pni) = — , 
(nlin)) = tripnj) = '-, 



Ani(7r) = 




In this particular example with Vm as the interaction, 
the above results actually is independent of r. Thus 
the measurement of total single particle state population 
fluctuations after the time reversed evolution provides a 
strong indication for the coherence and correlation of the 
created entangled states. 

VI. LIMITATIONS AND ADVANTAGES 

There are several limitations to our model that needs 
to be discussed carefully. First, the molecular coupling 
needs to dominate over the atomic coupling; i.e., we need 
to achieve a molecular coupling 

» - (75) 
Aa 

where Am /a s^re detunings from the molecular and 
atomic intermediate states, respectively. We can esti- 
mate analytically the free-bound transition amplitude us- 
ing the interacting ground state wave function from , 
and excited state wave function from the bound state so- 
lutions in a — Cs/r'^ potential modified by an inner repul- 
sive core 0. Most hkely, however, the deepest molecu- 
lar bound state with significant amplitude should be cho- 
sen in order to maximize the detuning from the atomic 
transition. In addition, we should have Am 3> 7a/ to 
minimize spontaneous emission which would lead to de- 
coherence and loss of atoms. 

A clear advantage of our scheme is the signature 
for diffraction via a molecular intermediate state. The 
diffracted atoms would move at half the speed of atoms 
that have undergone atomic Bragg diffraction. Even as 
the number of atoms is increased, our scheme achieves 
the same level of squeezing and the same high value of 
overlap with the massive GHZ state. 

Our scheme is based on an engineered interaction that 
can be turned on and off like the scheme of Sorensen and 
Molmer (Eq. (p|)). In particular, our scheme works for 
non-interacting atoms which can decrease the noise due 
to atom- atom interactions in a U(l) symmetry breaking 
condensate state. We have also provided explicit calcula- 
tions that demonstrate that the time reversed dynamics 
can be used to verify the coherence of the intended max- 
imally entangled states. 

Among all models discussed, there exists an even/odd 
number problem: The form of the final state depends on 
whether the total number of condensed atoms is even or 
odd number. If, however, the condensate resembles a co- 
herent superposition of number states (grand canonical 
ensembles), we should average the results over the num- 
ber distribution p{N) as well. This will be explored in a 
future study. 
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